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ABSTRACT 

It has been suggested for quite a long time that galaxy mergers trigger activities of supermassive black 
holes (SMBHs) on the grounds of imaging observations of individual galaxies. To quantitatively examine this 
hypothesis, we calculate quasar luminosity functions (LFs) by manipulating the observed galaxy LFs (z < 2) 
and theoretical merger rates from semi-analytical formulations. We find that the model reproduces the observed 
quasar LFs provided that the mass ratio (q) of the secondary galaxy to the newly formed one changes with 
cosmic time. The results show that the fraction of major mergers decreases from / ma j ~ 0.2 at z ~ 2 to / ma j — » 
at z ~ 0. As a consequence, the newly formed SMBHs from major mergers at z ~ 2 may acquire a maximal 
spin due to the orbital angular momentum of the merging holes. Subsequently, random accretion led by minor 
mergers rapidly drives the SMBHs to spin down. Such an evolutionary trend of the SMBH spins is consistent 
with that radiative efficiency of accreting SMBHs strongly declines with cosmic time, reported by Wang et al. 
(2009). This suggests that minor mergers are important in triggering activities of SMBHs at low redshift while 
major mergers may dominate at high redshift. 

Subject headings: galaxies: evolution — galaxies: high-redshift — quasars: general 
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1. INTRODUCTION 

In the hierarchical model of galaxy formation, quasar activ- 
ities are related to the mergers of galaxies (e.g., Toomre 1977; 
Hernquist 1989; Barnes & Hernquist 1992; Hopkins et al. 
2008). During the merging processes, tidal torque drives gas 
inflows into the inner region of galaxies, forming starbursts 
and igniting the central SMBHs (e.g., Hernquist 1989; Barnes 
& Hernquist 1996; Springel et al. 2005). Once quasar activi- 
ties are triggered, the mutual interaction between the SMBHs 
and galaxies is unavoidable accompanied by feedback of ac- 
tive galactic nuclei (AGNs) (e.g., outflows or radiation; see Di 
Matteo et al. 2005; Croton et al. 2006; Bower et al. 2006). 
This is believed to regulate coevolution of SMBHs and their 
host spheroids, inferred from the well-established correlation 
between SMBH mass and spheroid mass or dispersion veloc- 
ity (Magorrian et al. 1998; Gebhardt et al. 2000; Ferrarese 
& Merritt 2000; Marconi & Hunt 2003; Hiiring & Rix 2004), 
however, the details remain open. 

Based on the long-standing "merger hypothesis" (Toomre 
1977), quasar LFs have been extensively modeled within the 
framework of hierarchical structure formation (Roos 1985; 
Carlberg 1990; Haehnelt & Rees 1993; Haiman & Loeb 1998; 
Wyithe & Loeb 2002, 2003; Volonteri et al. 2003; Hopkins et 
al. 2008; Shen 2009). A brief summary of hierarchical mod- 
els from the published literature is given in Table 1 . We note 
that most models employ major mergers. In these models, 
the scaling relations between black hole mass and dark mat- 
ter halo mass are generally employed to combine the merg- 
ers of galaxies with dark matter halos. The dark matter ha- 
los follow the extended Press-Schechter (EPS) theory (e.g., 
Lacey & Cole 1993). After a merger event, a quasar is ig- 
nited and shines at the Eddington luminosity with a universal 
light curve. The rapid decreases as a main feature of quasar 
LFs can be well produced by the decrease in the merger rates 
and the availability of cold gas to fuel the SMBHs after the 



peak activities of SMBHs. Although these models quantita- 
tively describe the shape of quasar LFs, it is still insufficient 
to understand the cosmological evolution of quasars in light 
of merger hypothesis. 

There is increasing evidence showing that major mergers 
cannot constitute the only triggering mechanism of quasar 
activities from observations based on the studies of absorp- 
tion properties of AGNs (e.g., Hasinger 2008) or the galaxy 
morphologies (e.g., Georgakakis et al. 2009; Reichard et al. 
2009), as well as from theoretical work by modeling LFs of 
low-level AGNs (e.g., Hopkins & Hernquist 2006; Marulli et 
al. 2007). It is evident that minor mergers between galax- 
ies and smaller satellites occur much more frequently (Barnes 
& Hernquist 1992; Mihos & Hernquist 1994). Many pre- 
vious works have investigated the role of minor mergers in 
triggering quasar activities (Gaskell 1985; Hernquist 1989; 
Hernquist & Mihos 1995; De Robertis et al. 1998; Taniguchi 
& Wada 1996; Taniguchi 1999; Corbin 2000; Chatzichristou 
2002; Kendall et al. 2003; Hopkins & Hernquist 2009a). As- 
suming that growth of SMBH occurs via episodic accretion in 
light of Soltan's argument (Soltan 1982), Wang et al. (2009, 
hereafter W09) derived an ^-equation to describe cosmolog- 
ical evolution of radiative efficiency of SMBHs. It is found 
that the radiative efficiency does exhibit a strong cosmologi- 
cal evolution: it decreases from ij ~ 0.3 at z ~ 2 ! to r\ ~ 0.03 
at z ~ 0. This can be interpreted by the mounting importance 
of minor mergers towards low redshifts. At high redshift, ma- 
jor mergers produce fast rotating SMBHs due to conservation 
of (orbital) angular momentum of the merging holes. Sub- 
sequently random accretion, which may be related to minor 
mergers (Kendall et al. 2003; King & Pringle 2007), acts to 
spin them down efficiently. 

In this paper, we extend previous model (Wyithe & Loeb 
2002) and construct quasar LFs from the observed galaxy LFs 
and theoretical merger rates. Our goal is to investigate the 
redshift-dependent merger history and the resulting spin evo- 
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TABLE 1 
Model list 



Reference 


Model 


Merger type a 




Quasar light curve 


M, 


-Mh relation 






^-dependent 


Mil -dependent 


R85 


Star disruption + galaxy merger rate 


Yes 








C90 


Halo merger rate 


Major 


Gau 






HR93 


Halo formation rate 




Exp 


V 


X 


HL98 


Halo formation rate 




Exp 


X 


X 


WL02 


Halo mass function + halo merger rate 


Yes 


Step 


V 


v 7 


WL03 


Halo mass function + halo merger rate 


Major 


Step 


v 7 




VHM03 


Merger tree 


Major 


Step 


X 


< 


HHCK08 


Halo mass function + galaxy merger rate 


Major 


Pow 


v' 




S09 


JV-body merger rate 


Major 


Exp + Pow 


v' 




This work 


Galaxy LFs + halo merger rate 


j-dependent 


Step 


V 





REFERENCES. — R85: Roos (1985); C90: Carlberg (1990); HR93: Haehnelt & Rees (1993); HL98: Haiman & Loeb (1998); WL02: Wyithe & Loeb (2002); 
WL03: Wyithe & Loeb (2003); VHM03: Volonteri et al. (2003); HHCK08: Hopkins et al. (2008); S09: Shen (2009). 
a "Yes" means no specifying merger type. 

b "Exp" means exponential, "Pow" means pow law, "Step" means step function, and "Gau" means the Gaussian distribution of quasar lifetime.. 
c " v / " means dependent and "x" means independent. 



lution of SMBHs. We describe our model in Section 2. In 
Section 3, we present comparisons of the modeled quasar LFs 
with the observational data, followed by discussions on the 
merger history and accretion history, and the spin evolution of 
SMBHs. The conclusions are drawn in Section 4. We use cos- 
mological constants VL m = 0.3, f?A = 0.7, the Hubble constant 
Hq = 70 km s" 1 Mpc" 1 and a% = 0.87 throughout the paper. 

2. MERGERS AND BLACK HOLE ACTIVITIES 

2.1. Merger rates of galaxies 

It is generally assumed that galaxies reside in the central 
region of dark matter halos, and consequently galaxy mergers 
follow dark matter halo mergers. In our calculations, this is 
implemented through the relationship between the masses of 
central SMBHs in galaxies and their host dark matter halos. 
First, we present merger rates of dark matter halos based on 
the EPS theory in this section. The probability for a halo of 
mass Mi at time t merging with another halo of mass AM is 
given by (Lacey & Cole 1993) 
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where M2 = Mi + AM is the mass of newly formed halo, 01 = 
<t(Mi), and cr 2 = o{Mi). The expression of a is given by 



a\M) = 
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where j(x) = (sinx-xcosx)/* 2 , p m = 3// 2 0,„/87rG, and M = 
4-Trp m R 3 /3 is the mass within a comoving sphere of radius R, 
G is the gravity constant, H is the Hubble constant, and Vl m is 
the fraction of baryons. The primordial power-law spectrum 
is 

■PokT 2 (k), (3) 

-1/4 



P(k): 



where the transfer function is T(k) = yl + J2i=i c i1' 

(c s qy 1 \n(l + c 5 q), ci = 3.89, c 2 = 259.21, c 3 = 162.77, c 4 = 
2027. 17, C5 = 2.34, po is determined by observations of a% at 



8/z 1 Mpc and^ = A:/(r2 m /j 2 Mpc '), where h is the Hubble con- 
stant in a unit of 100 km s" 1 Mpc -1 . 5 c (t) is expressed by 

2/3 /fn\ 2 / 3 



: S c0 [ J 



91 7 



(4) 



S c (t) ■ 

where to is a referenced epoch. 

2.2. Quasar luminosity functions 

In modeling quasar LFs, we extend the method of Wyithe & 
Loeb (2002) by employing the realistic galaxy LFs instead of 
the halo mass function and considering the redshift-dependent 
merger history. To derive quasar LFs from known galaxy LFs 
through mergers, we assume that (1) each galaxy hosts an 
SMBH; (2) all mergers trigger SMBH activities; (3) a merger 
simultaneously triggers an SMBH activity. We neglect the 
time delay between a galaxy merger and igniting of the cen- 
tral SMBH. 

We use the scaling relation between the black hole mass 
and the dispersion velocity M. oc er 4 (Tremaine et al. 2002) to 
connect SMBHs with dark matter halos. The dispersion ve- 
locity correlates with the Keplerian velocity of the dark matter 
halo via v c oc cr' 3 with /3 ~ 0.84 (Ferrarese 2002), where v c is 
determined by Mh (Barkana & Loeb 2001) 



v c = 159.4M, 
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where Mh.h = Mh/10 hM® is the dark matter halo mass, 
A, = 187r 2 + 82t/-39<i 2 is the final over-density relative to 
the critical density at redshift z, d = fl m . z - 1 and VL m = 
[Q m (l +z) 3 ] I [O m (l +z) 3 + fi A ] ■ The ratio of SMBH mass to 
halo mass is then computed as (Wyithe & Loeb 2002) 



— = eoM£ u 
Mh 



A c 



1/6 
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where 7 = 4//3 = 4.76 and eo is a free parameter. We note 
that e is weakly dependent on M H , but sensitive to redshift. 
The nonlinear relation between SMBH mass and halo mass 
with 7 > 3 is due to the significant baryon accretion during 
SMBH activities visible to observers (Yu & Tremaine 2002). 
It is worth pointing out that the potential evolution of 7 is 
poorly understood (e.g., Shields et al. 2003; Jahnke et al. 
2009; Shields et al. 2006; Alexander et al. 2008; Coppin et 
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TABLE 2 

The free parameters in the models 



Parameter 


Given by galaxy LFs of a 


Implications 




llbert et al. (2005) 


Dahlen et al. (2005) 




log(eo) 


-3.3 ±0.7 


-5.0 ±0.6 


The initial ratio of halo to SMBH mass 


log(?dc,o/yr) 


7.54 ±0.04 


7.30 ±0.04 


Characterized timescale of a single episode 


A 


0.10±0.01 


0.20 ±0.01 


The mean Eddington ratio of quasars 


q(.z) 






AM. /M. as a function of z 



NOTE. — The xj,j n /d.o.f. of fittings are 78/65 and 121/65 for llbert et al. and Dahlen et al., respectively. 
a Errors are calculated by assuming A\t /'J = 0.1. 



al. 2008). In our model, we take 7 = 4.76 and neglect the 
potential effects of its evolution. The modeled quasar LFs are 
found to be insensitive to 7. 

Once a quasar is triggered, it is assumed to shine with a uni- 
versal B-band light curve as (see also Haiman & Loeb 1998) 

L B (t)=M.f(t) = eM u f(t), (M H >M H . mm ) (7) 

where /(f) is a function of time and M H m ; n is the minimum 
mass of the halo, inside which an SMBH can form. Here we 
simply use a step function for /(?) (Wyithe & Loeb 2002), 



AWAMh 
J C B M. \ M H ■ 



(8) 



where L Ed d = 1.38 x 10 (M./M Q ) is the Eddington luminos- 
ity, A is the Eddington ratio, AMh is the mass of the secondary 
halo, O(f) is the Heaviside step function, Cb is the bolo- 
metric correction factor defined as the bolometric luminosity 
^Boi = CbLb (Marconi et al. 2004). This light curve takes into 
account that the gas deposited from the satellite galaxy runs 
out in a time proportional to AMh /Mh Q (Wyithe & Loeb 



1 Basically, the gas from the satellite galaxy is 
AM. 



AAfgas AM H . . AM H . . 

Mh <x — — — M., 



AMh Mh ** Mh 
where the ratios e = M. /Mh and AMg as /AMn are assumed as constants in 



2002). It turns out that spectral energy distributions (SEDs) of 
AGNs have no appreciable evolution with cosmic time (e.g., 
Shemmer et al. 2008). This would guarantee that Eddington 
ratio (A) may also be independent of redshift, according to its 
correlation with the SEDs (Wang et al. 2004). A similar con- 
clusion has been drawn through an analysis of a large quasar 
sample from the Sloan Digital Sky Survey (SDSS; Shen et 
al. 2008) and the AGN and Galaxy Evolution Survey (AGES; 
Kollmeier et al. 2006), respectively. Therefore, we assume 
that the triggered SMBHs shine at a fixed Eddington ratio. 
The value of A can be obtained from fittings of quasar LFs. 

With these assumptions, the modeled quasar LFs can be 
written as (Wyithe & Loeb 2002) 

r qM. 



dz' AM. S[L B -M.f(t z -t')] 



d 2 N 

u JV merg 




dn. 




dAM.df 




dM. 





— dAM., 



(9) 



where t z and t 1 are the cosmic time at redshift z and z', S(x) is 
the Dirac delta function, |dn./dM. \ M , is the number density 

of black holes, |(d 2 Af me i- g ) / (dAM.df) I is the merger rate, 

a zeroth-order approximation. As accretion of the SMBH is limited at the 
Eddington rate, the gas will run out via accretion in a time proportional to 
(AMh/Mh). 
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FIG. 2. — (a) Cosmological changes of the mass ratio (q), and (b) the frac- 
tion of major mergers to total (/,-•) as a function of redshift, indicating the 
increasing role of major mergers with redshift. Solid and dashed lines (with 
symbols) are for Ilbert et al. (2005) and Dahlen et al. (2005), respectively. 
Error bars are calculated by assuming = 0.1. 

M' m =M % - AM. and the mass ratio 



q(z)- 



AM. M.-M'. 



M. 



M. 



(10) 



Note that a major merger is defined as 0.5 > q > 0.25 so 
that the mass ratio of the two merging galaxies is larger than 
1:3. The upper limit forces the secondary galaxy to always 
be smaller than the primary one. We treat q as a function of 
redshift and constrain it from the fittings of quasar LFs. 
Integrating overM., we can find 



dz' 



M.f(t z -t') 7 



rqM. 


d 2 N 

u 1 v merg 


Jo 


dAM.df 



dn. 



dM. 



dAM.. 



(11) 



Quasars in the cosmic time range (f z -fg, f,) with a lifetime 
tQ = (AMe/Me^dcO contribute to the LFs at cosmic time f-. 
Since the lifetime is much smaller than the cosmic time f- ~ 
H(z)~ l at the redshift z < 2, we can then integrate over z' by 
replacing dt 1 with (AMn/Mu)tdc.o- We employ the galaxy 
LFs, defined as $(L',z), to trace the SMBH mass function 
dn./dM., and use the ratio (M./Mh) to calculate the merger 
rate of SMBHs from that of dark matter halos, 



d 2 /V 

u JV merg 




3 


d 2 N 

u 1 * merg 


dAM.df 




7e 


dAM H df 



7e 



(12) 



where M' H = M H - AM H . Finally we get the quasar LFs 



dL' T H • / AM H a 
nL B ,z) = t dc , Q — / K$(L',z)— -dAM H , (13) 

dl^B Jo.05M H m h 
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FIG. 3 . — (a) Duty cycle time of black hole activity as a function of quasar 
luminosity at the redshift bins corresponding to Figure 1 from bottom to top. 
(b) Luminosity-averaged duty cycle time with redshift. Error bars are the up- 
per and lower limits of the duty cycle time in (a). Solid and dashed lines (with 
symbols) are for Ilbert et al. (2005) and Dahlen et al. (2005), respectively. 



where the factor (dL'/dL^) comes from the different lumi- 
nosity bands used in the galaxy and quasar LFs, and L' is 
the luminosity of galaxies hosted in the halos with a mass 
Mh - AMh- We calculate L' using the relation between the 
/^-magnitude of normal galaxies and the mass of black holes 
from McLure & Dunlop (2002; see the appendix). Note that 
the low limit in Equation (13) is set so that the secondary 
galaxy should be massive enough to provide significant gas 
for the mass assembly of the newly formed black hole. 

With the merger rates, we have a duty cycle time Q of the 
quasars during a Hubble time H~ 



[ qMn AM H , • . 
tdc(L B , z) = fdc.o / -77^H: x TZd AM H , 

0.05M H M H 



(14) 



where H, is the Hubble constant at redshift z- We also define 
a fraction of major mergers to the total as 



/majfcM H )=-7 



AUj(z,M H ) 



where 



AUz,M H ) 



(15) 



AUj(z,M H )= / TZ(z)dAM H , (16) 

J0.25Mh 

2 Note that here "duty cycle time" refers to the active time of the quasars 
in a timescale ffr 1 . 
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and 

rqM H 

NUz,Mn)= / 7t(z)dAM H . (17) 

J0.05M H 

According to the definition of a major merger, when q{z) < 
0.25, we have / ma j = 0. We stress here that the parameter 
/maj(z) conveniently traces the gradual evolution from major 
to minor mergers, showing a merger history. 

As a summary, Table 2 gives a list of the free parameters 
used in our model. The parameter q determines both the shape 
and the number density of quasars. The amplitude of quasar 
LFs is linearly proportional to f ( j c ,o ! which is easily justified 
from Equation (13). The other parameters (eo and A) can be 
independently quantified by theoretical or observational argu- 
ments. Fitting quasar LFs would provide more precise values 
of these parameters deliberately. 

3. COMPARISON WITH OBSERVATIONS 

3.1. Modeling quasar LFs 

We use two sets of galaxy LFs well-determined by Hubble 
Space Telescope and the Great Observatories Origins Deep 
Survey (GOODS; Dahlen et al. 2005), and VIMOS-VLT 
Deep Survey (VVDS; Ilbert et al. 2005), respectively. Since 
both sets of galaxy LFs are only available for z < 2.0, the 
modeled quasar LFs are then also limited. Galaxy LFs are 
generally described by, 

c^) = c^£) %(-£), (18) 

where $» is a constant, a is an index and L» is the cutoff 
luminosity of galaxies. We have to keep in mind that the 
galaxy LFs are based on photometric redshifts. We use quasar 
LFs determined by the combined data of the Two-Degree 
Field survey and the SDSS, which to date are the best con- 
strained LFs with unprecedented precision (10,637 quasars) 
and the dynamic range (g-band flux limit of 21.85; Croom et 
al. 2009). 

Quasar LFs are characterized by double power-laws with a 
break luminosity (e.g., Hopkins et al. 2007). In our model, 
the break luminosity is mainly determined by the cutoff lu- 
minosity of galaxies (L»). The indices of the double power- 
laws of quasar LFs are jointly determined by a and the merger 
rates. In modeling quasar LFs, the free parameters (A, eo and 
f<k,o) degenerate weakly. However, their typical values can be 
quantified directly from observations. We fixA~ 0.15 (Shen 
et al. 2008), f dc ,o ~ 3.0 x 10 7 yrs (Martini & Weinberg 2001), 
and eo ~ 10~ 5 (Ferrarese 2002). We calculate a grid of quasar 
LFs in a range of A, ?dc,o and eo around their typical values and 
then obtain the fitting values accordingly. 

The best fittings with minimal \ 2 are given in Table 2 and 
shown in Figure 1 . We find that the global features of quasar 
LFs can be well produced by Equation (13). The results mod- 
eled by galaxy LFs from Ilbert et al. (2005) and Dahlen et 
al. (2005) are in reasonable agreement. To match the obser- 
vational data, we find that q must be a function of redshift for 
both sets of galaxy LFs. We get a dependence of the mass 
ratios q with respect to redshift as shown in Figure 2(a). It is 
found that q exhibits a rapid decline with redshift. A more ex- 
plicit parameter f ma ; clearly shows significant evolution with 
redshift in Figure 2(b): it decreases from / ma ; ~ 0.2 at z <~ 2 to 
/maj — ► at z ~ 0. This suggests a merger history with a trend 
from major to minor since z ~ 2. Major merger rates mea- 
sured from the Hubble Ultra Deep Field get a peak at z <~ 1-2 



(Ryan et al. 2008), further supporting the evolution of / ma j 
obtained here. 

The duty cycle time is obtained from Equation (14), and 
the results are illustrated in Figure 3(a). We find that fdc is 
strongly dependent on redshift, but insensitive to the quasar 
luminosity, namely, the SMBH mass. The typical duty cycle 
time is > 10 7 years, indicating that the SMBHs are undergo- 
ing many episodes (Martini 2004; Wang et al. 2006; 2008; 
2009). The lifetime of quasars is longer at higher redshift. 
This should be due to the increase of merger rates and the 
mounting importance of major mergers with increasing red- 
shift. Major mergers provide more abundant gas, and thus 
drive SMBH activities to live longer than minor mergers. We 
also calculate the luminosity-averaged duty cycle time in Fig- 
ure 3(b), and simply fit its relation to redshift as 

f dc (z) = 2xl0 6 (l+z) 2 ' 5 yr. (19) 

The evolution of the duty cycle time should be independently 
examined by future observational data, for example from the 
clustering measurements of quasars (for a review see Martini 
2004). In addition, a more realistic model of the quasar life- 
time and light curve (Hopkins et al. 2005; Hopkins & Hern- 
quist 2009b) is worth considering in a future paper. 

For the bin of 0.4 < z < 0.68, we find the modeled LFs 
under and overestimate the observed number of quasars at 
faint and bright tails, respectively. Bonoli et al. (2009) 
also find similar behaviors from semi-analytically modeling 
quasar LFs coupled with numerical simulations. The over- 
estimate for the bright tail is because most of the galaxies are 
gas-poor (dry mergers, see Lin et al. 2008) and some fraction 
of (low-level) quasar activities may be stochastical accretion- 
driven (Hopkins & Hernquist 2006). This naturally drives an 
underestimate of quasars for the faint tail regarding passive 
evolution of galaxies after z ~ 1 (Faber et al. 2007). 

3.2. Mergers and accretion history 

Major mergers are undoubtedly taking place in very lumi- 
nous submillimeter galaxies at z ~ 2-3 (Tacconi et al. 2006; 
2008). Detailed studies from Spectroscopic Imaging survey 
in the near-infrared with SINFONI (SINS) of high redshift 
galaxies (z <~ 1-3) show that about one third of the galax- 
ies are rotation-dominated rings/disks undergoing star forma- 
tions, another one third are compact and velocity dispersion- 
dominated objects, and the remaining are clearly interact- 
ing/merging systems (Genzel et al. 2008; Forster Schreiber 
et al. 2009). It appears that rapid and continuous gas accre- 
tion via "cold flows" and/or minor mergers are playing an im- 
portant role in driving star formation and the assembly of star 
forming galaxies (Genzel et al. 2008). Based on morpholo- 
gies of X-ray selected AGNs with redshifts 0.5 < z < 1.3, 
Georgakakis et al. (2009) suggested that AGNs in disk galax- 
ies are likely fueled by minor mergers. More interestingly, 
Robaina et al. (2009) reported that only less than 10% of star 
formation in massive galaxies for z ~ 0.6 is triggered by major 
interaction. These pieces of independent evidence are consis- 
tent with the change of the major merger fraction as shown in 
Figure 2b. 

The importance of minor mergers in triggering SMBH ac- 
tivities has been realized by many authors (e.g., De Robertis et 
al. 1998; Hernquist & Mihos 1995; Taniguchi 1999; Kendall 
et al. 2003; Hopkins & Hernquist 2009a). The dynamical 
studies of minor mergers find that the merging satellites re- 
tain memory of the initial random orientations when arriving 
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at the nucleus of the primary galaxy (Kendall et al. 2003). We 
therefore may infer that during a succession of minor merg- 
ers, gas fueling to the central SMBH will carry on in a man- 
ner with stochastical orientations relative to the host galactic 
plane. This affords a natural explanation of random accretion 
proposed by King & Pringle (2006; 2007), as one fueling sce- 
nario to AGNs. There are several independent indications for 
minor mergers as a triggering mechanism of random accre- 
tion. 

First, radio jets in a sample of nearby Seyfert galaxies 
(z < 0.03) have a random orientation with respect to their host 
galactic plane (Kinney et al. 2000). Presumably, the jet direc- 
tion aligns with the axis of the accretion disk, therefore the fu- 
eling gas to the SMBHs could not be supplied from their host 
galaxies directly (see a brief review of fueling mechanisms of 
SMBHs in Taniguchi 1999), but probably from minor mergers 
(King & Pringle 2007). 

Second, AGN types are independent of orientations of their 
host galaxies (Keel 1980; Dahari & De Robertis 1988; Munoz 
Marin et al. 2007). Additionally, the orientations of extended 
[O III] emission in some nearby Seyfert galaxies have no cor- 
relation with the major axis of the host galaxies (Schmitt et al. 
2003). This suggests that the dusty torus axis is not aligned 
with the host galaxy rotation axis. 

Third, W09 derived an //-equation to describe cosmologi- 
cal evolution of radiative efficiency of accretion into SMBHs, 
based on the argument that the SMBH growth is dominated by 
baryon accretion (Soltan 1982). Applying the equation into 
existing survey data of quasars and galaxies, they report that 
the radiative efficiency strongly evolves with redshift, with 
values of 77 ~ 0.3 at z ~ 2 and 77 ~ 0.03 at z ~ 0. As de- 
scribed below, episodic random accretion definitely leads to 
a decrease of spin for SMBHs. Presumably, with the link be- 
tween radiative efficiency and SMBH spin (Thorne 1974), this 
trend clearly supports the cosmological evolution of mergers 
from major at high-z to minor at low-z found from our model 
of quasars LFs. 

With the above comprehensive indications and the results 
obtained in our model of quasar LFs, it quite appears that 
minor mergers may mainly be driving SMBH activities af- 
ter z ~ 2.0, and random accretion driven by minor mergers 
supplies shining power of the quasars. 

3.3. Spin evolution 

During the active phases of a quasar, the SMBH accretes a 
significant amount of gas to buildup its mass (the Soltan's ar- 
gument; see also Yu & Tremaine 2002; Volonteri et al. 2003). 
The previous studies using hierarchical merger tree found that 
mass accretion also plays a dominant role in determining spin 
evolution of black holes compared with coalescence (Volon- 
teri et al. 2005; Berti & Volonteri 2008). However, the de- 
tailed spin distribution of black holes depends on the specific 
accretion scenarios: prolonged or random accretion. In pro- 
longed accretion, the material keeps constant sign of angu- 
lar momentum. In this case, an initially nonrotating black 
hole that accretes about one-third of its mass will end up very 
fast spinning (Thorne 1974). By contrast, random accretion 
rapidly spins down black holes (Berti & Volonteri 2008). On 
the other hand, in a major coalescence some orbital angular 
momentum of a black hole binary may be converted to the 
spin of the newly formed black hole. This most likely yields 
a rapidly rotating black hole. 

As mentioned above, random accretion can be naturally 
interpreted through stochastically-oriented minor mergers 



(Kendall et al. 2003). We have shown that minor mergers 
become important in driving SMBH activities toward low red- 
shifts. Then we can estimate spin evolution of SMBHs due to 
random accretion induced by minor mergers. For an approx- 
imation, we apply the equation a = ao(M,x)/M,) 2A (Equation 
(26) in Hughes & Blandford 2003, for random captures of 
minor black holes originally). Here a = JcjGM\ is the spin 
parameter of a black hole with a mass M, and angular mo- 
mentum /, c is the speed of light, aq and M. are the initial 
spin and mass of the hole, respectively. Under the assump- 
tion that each merger triggers a quasar activity, the number of 
merger events can be estimated through the ratio of the duty 
cycle time to the average lifetime of a single activity 



tic 



AfQ ' 



(20) 



where A?q is the average lifetime of a single activity, tdc is 
given by Equation (19). The lifetime of a single activity can 
be inferred from the proximity effect, AfQ ~ 1 Myr (Kirkman 
& Tytler 2008). For z ~ 2, we apply f dc « 3 x 10 7 yr and then 
obtain the order of magnitude estimate n <~ 10. To satisfy the 
ratio of SMBH mass and halo mass in Equation (6), the mass 
growth of SMBHs due to accretion during each merger is (see 
also Volonteri et al. 2003) 



Am = ^AM H = ~(0M., 



(21) 



where (£) is the average mass ratio of two merging galax- 
ies/halos. The total mass growth during a Hubble time is thus 



AM. nAm 7 

= = —n(t) ~ 2, 

Mm M. 3 vs/ 



(22) 



where 7 = 4.76 and we set (£) ~ 0.1 for minor mergers. Con- 
sequently, the spin of the black hole drops off by a factor of 
3 2 4 ~ 10 from z ~ 2 to z ~ 0. This is generally consistent with 
the results found by Berti & Volonteri (2008) that random ac- 
cretion predicts the spin distribution with a < 0.1. 

We may suggest a picture that at high redshift z ~ 2, the 
newly formed SMBHs from major mergers acquire a rapid 
spin due to the orbital angular momentum of the merging 
holes, and subsequently random accretion led by minor merg- 
ers may rapidly spin them down. It is interesting to note that 
King et al. (2008) proposed a scenario that self-gravity of a 
accretion disk can give rise to a long series of episodic ac- 
tivities during a major merger event. We neglect this effect 
in our model. A future work combining spin evolution with 
semi-analytical theory will provide new insights into the co- 
evolution among SMBHs, galaxies and dark matter halos. 

4. CONCLUSIONS 

With the assumption that black holes are triggered through 
galaxy mergers, we successfully reproduce the observed 
quasar LFs by manipulating the galaxy LFs and merger rates 
from semi-analytic theory from z ~ to z ~ 2.0. This sup- 
ports that mergers are driving black hole activities with many 
episodes. We find that minor mergers become important in 
triggering quasar activities towards low redshifts. As a natu- 
ral consequence of minor mergers, the random accretion may 
drive growth and spin evolution of SMBHs at low redshift. 
The estimated spin declines rapidly with decreasing redshift, 
in agreement with the evolutionary trend of radiative effi- 
ciency for SMBHs. 
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TABLE 3 

The relationship log(M./M ) = a+bM R from published literature 



Reference 


Sample" 


Band b 


a 


b 


Scatter 


Note 


MD02 


72S+18E 


R 


-2.96 ±0.48 


-0.50 ±0.02 


0.39 




BFFF03 


20E 


R 


-3.00 ±1.35 


-0.50 ±0.06 


0.39 




MH03 


12E+10S 


8 


-2.54 ±1.10 


-0.50 ±0.05 


0.48 




EGC04 C 


8E+5S 


R 


+2.70 


-0.25 


0.35 




G07 


12E+10S 


R 


+1.40 ±0.33 


-0.38 ±0.04 


0.33 


sample bias? 


L07 


100E+55S+64B 


V 


-3.70 ±0.77 


-0.55 ±0.03 




K08 


45S 


R 


-2.87 ±0.04 


-0.50 


0.41 


fix/3 


BPPV09 


19E 


V 


-1.67 ±1.58 


-0.47 ±0.08 


0.64 


fitting routine? 


BPPV09 


26A 


V 


+1.22±0.75 


-0.32 ±0.04 


0.44 


fitting routine? 



REFERENCES. — MD02: McLure & Dunlop (2002); BFFF03: Bettoni et al. (2003); MH03: Marconi & Hunt (2003); EGC04: 
Erwin et al. (2004); G07: Graham 2007; L07: Lauer et al. (2007); K08: Kim et al. (2008); BPPV09: Bentz et al. (2009). 
a "E" is for elliptical galaxies, "S" is for either SO galaxies or Spiral galaxies, "A" is for AGNs, and "B" is for bright cluster 
galaxies. 

b B—R = 1.52 and V—R = 0.61 from Fukugita et al. (1995) are used to convert B- and V-band magnitude to S-band magnitude. 
The standard relation logLy /Lq = 0.4(-My +4.83) is employed to calculate My from L v (Bentz et al. 2009) and M b .q = 5.47 is 
used (Cox 2000). 

c We artificially extract the data from the Figure (1.3) in EGC04. 



We thank the anonymous referee for the helpful com- 
ments and suggestions that greatly improved the paper. YRL 
acknowledges H.-B. Yuan and W.-W. Zuo for reading the 
manuscript and I. A. McNabb for the help in improving the 
English presentation. We also appreciate the stimulating dis- 
cussions among the members of IHEP AGN group. The re- 
search is supported by NSFC-10733010 and 10821061, CAS- 
KJCX2-YW-T03 and 973 project (2009CB824800). 

APPENDIX: INFLUENCES OF M.-L BUL RELATIONSHIP ON QUASAR 

LFS 

The relationship between SMBH mass and bulge luminos- 
ity (M.-Lb u i) is an important ingredient in modeling quasar 
LFs from the galaxy LFs. Many previous works have focused 
on this relationship (see Table 3 for a brief summary). Unfor- 
tunately, not all of them reach the consistent conclusions. At a 
first glance, the quiescent galaxies have a steeper slope of M.- 
Lbui than the AGNs. This may be due to neglecting the effect 
of radiation pressure, which causes underestimations of black 
hole mass in AGNs as proposed by Marconi et al. (2008, see, 
however, Netzer 2009). Meanwhile, Graham (2007) recently 
made an attempt to compensate for the analysis techniques 
employed by previous groups (McLure & Dunlop 2002; Mar- 
coni & Hunt 2003; Erwin et al. 2004). They took into ac- 
count a number of issues, e.g. the removal of a dependency 
on the Hubble constant and a correction for the dust atten- 
uation in bulges and disk galaxies. After careful revisions 
and adjustments, they found the differences among previous 
groups could be eliminated and they reproduce a somewhat 
shallower slope. However, they did not exclude the potential 
sample bias, which may lead to a similar relation as McLure 
& Dunlop (2002) [see Equation (A7) in Graham 2007]. 

To sum up, there appear to be many unclear aspects in the 
M,-L\, u \ relationship that deserve further examinations with 
improved observational data. We adopt McLure & Dunlop's 
relation in the present work. To show the influences of M.- 
Lbui relationship on the quasar LFs, we simply perform an ap- 
proximate estimate. For AMh/Mh C 1, we have the merger 



rate (Lacey & Cole 1993) 



/AMh 



-3/2 



(23) 



Using the M. -Lbui relation and the mass ratio between the 
black hole and the dark matter halo, we have 



M H = M./eoc 10 hM "/e, 
and the galaxy LFs in /?-band 

$(M«) oc l0 a4(1+a >( M «*- M «)exp Lio a4(M «- M «> 



(24) 



M h 



0.4(l+a)/6 



exp 



M h 



0.4/6' 



(25) 



where M H = l0 a+hM '< /e. Combining the above equations into 
Equation (13) and performing some mathematical expansions 
for AM H /M H <C 1, we have the quasar LFs approximately as 



* oc 



-1/2 



x exp 



Mi 



L* 



M H /M H 
1 - AMh/Mh 

0.4(l+a)/i 

exp 



M*/M H 
1 - AMh/Mh 

0.4/fc' 



0.4(l+a)/6 



Mpi 
M H 



dAM H 

0.4/6" 



0.4(l+a)/i> 



exp 



0.4/6' 



(26) 



We can find that the break luminosity L* of the quasar LFs is 
dependent on the cutoff luminosity M% of the galaxy LFs, and 
the slope 0.4(1 +a)/b is relative to the slope of the M.-Lbui 
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